Optimal rotations of deformable bodies 
and orbits in magnetic fields 
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Deformations can induce rotation with zero angular momentum where dissipation is a natural 
"cost function". This gives rise to an optimization problem of finding the most effective rotation 
with zero angular momentum. For certain plastic and viscous media in two dimensions the optimal 
path is the orbit of a charged particle on a surface of constant negative curvature with magnetic 
field whose total flux is half a quantum unit. 



Rotations with zero angular momentum are intriguing. 
The most celebrated phenomenon of this kind is the ro- 
tation of a falling cat. A mechanical model Q replacing 
the cat by two rigid bodies that can rotate relative to 
each other, has been extensively studied, see 0, 01 ancl 
references therein. Here we address rotations with zero 
angular momentum under linear deformations. Our mo- 
tivation comes from nano-mechanics: Imagine an elastic 
or plastic material with its own energy source, and ask 
what is the most efficient way of turning it through an 
appropriate sequence of autonomous deformations with- 
out external torque. 

Deformations can generate rotations because order 
matters: A cycle of deformations will, in general, result in 
a rotation. The limiting ratio between the rotation and 
the area of the controls (deformations), when the lat- 
ter tends to zero can be interpreted as curvature 0, . 
Consequently, small cycles are ineffective since a cycle of 
length s in the controls leads to a rotation of order e 2 . 
The search for optimal paths forces one to mind defor- 
mations that are not small. 

The problem we address has three parts. The first 
part is to determine the rotation for a given path of the 
controls. We solve this problem for general linear defor- 
mations. In two dimensions this leads to curvature on the 
space of the controls which is localized. If one thinks of 
the curvature as a magnetic field, then the total magnetic 
flux is half the quantum unit. The second part is to set 
up a model for the cost function which we choose to be 
a measure of dissipation. We focus on two settings, one 
where the dissipation is rate independent, as is the case 
in certain plastics, and the other where it is rate depen- 
dent as in liquids. Both cost functions lead to the same 
metric on the space of deformations. The third part is to 
pose and solve the problem of finding the path of mini- 
mal dissipation for a given rotation. In two dimensions 
and for either model of dissipation, the problem maps to 
finding the shortest path that starts at a given point and 
encloses a given amount magnetic flux. As we shall see, 
optimal paths tend to linger near the circle in the space of 
controls where the ratio of eigenvalues of the quadrupole 
moment of the body is (4> + \f$) 2 . <f> = (1 + is the 

golden ratio. 



Deformations generate rotations because angular mo- 
mentum is conserved ja]- Consider a collection of point 
masses m a at positions x a . Internal forces may deform 
the body, but there are no external forces. Suppose that 
the center of mass of the body is at rest at the origin 
and that the body has zero angular momentum. The to- 
tal angular momentum, Ly, must then stay zero for all 
times. 

A linear deformation is represented by a matrix M 
that sends x — > Mx. The i,j component of the angular 
momentum is 



Li 



TV (MQM* £ij), 



where Q is the quadrupole moment of the body 



(i) 



(2) 



and £ij, i < j are the generators of rotations in n 
dimensions, i.e. (x, iijy) — Xiy.j —XjUi. Since tij span the 
anti-symmetric matrices, the set of equations L,j = 
imply that the matrix {dM)QM t is symmetric. 

Two immediate consequences of this symmetry are: 

• Isotropic bodies: With Q = 1, and M close to the 
identity, the symmetry of (dM)QM t reduces to dM 
being symmetric: the linear transformation must 
be a strain Q. 

• Pointers: Pointers are bodies with large aspect ra- 
tios, such as needles and discs. In the limit of infi- 
nite aspect ratio, Q may be identified with a pro- 
jection where dimQ is the dimension of the pointer. 
With M near the identity, the symmetry of (dM)Q 
implies that (1 - Q)(dM)Q = (1 - Q)Q{dM) t = 0. 
Since Q does not acquire a component in the nor- 
mal direction, 1 — Q, under dM a pointer keeps its 
orientation. 

We now derive the fundamental relation between the 
response (rotation) and the controls (deformations). To 
this end we use the polar decomposition M — RS with 
R orthogonal and S positive. Assuming S positive is a 
choice of a gauge which makes the representation unique 
with 5" = \jM l M. The symmetry of (dM)QM t gives 

{A, SQS} = SQ(dS) - (dS)QS, A = R^dR. (3) 
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Eq. © determines the differential rotation, A, in terms of 
the variation of the controls dS. The symbol d stresses 
that the differential will not, in general, integrate to a 
function on the space of deformations. Geometrically, 
the differential rotation is the connection 1-form which 
fixes a notion of parallel transport. 

Eq. (|3J) can be interpreted in terms of a variational 
principle: The motion is such that the kinetic energy is 
minimal for a given deformation. To see this, let M = 1 
and M = R + S with R antisymmetric (i.e. a rotation) 
and S symmetric (i.e. a strain). The kinetic energy is 
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f 



E = -Tr (MQM*) = - Tr ((R + S)Q(-R + S)) . (4) 
Minimizing with respect to R gives 

= SE = ^Tr (5R( - {Q, R} + [Q, S}) . (5) 

The trace is of a product of antisymmetric matrices, and 
its vanishing for an arbitrary antisymmetric SR implies 
{R, Q} = [Q, S] which is Eq. © for M = 1. 

One readily sees that if A is the solution of Eq. © 
given S and Q, then it is also a solution for XS and Q 
for A a scalar valued function. Hence scaling does not 
drive rotations and we may restrict ourselves to volume 
(or area) preserving deformations with det S = 1 . 

Since any Q is obtainable by a linear deformation of 
the identity, we may assume without loss of generality 
that Q = 1. Eq. reduces to 



{A, S 2 } = [S, dS] 



(6) 



Eq. @ is conveniently solved in a basis where S is diag- 
onal. Let Sj denote the eigenvalues of S then 



•Ail 9 , 9 (dS^i 

sf + s] 



(7) 



The curvature F is more interesting than the connection. 
It is defined by F = dA + A A A. Calculation gives 

^ (s 2 + s])(s] + s 2 k )(s 2 + si) 

The situation is particularly simple in two dimensions. 
We use the Pauli matrices a x ,a z , ia y and the identity as 
basis for the real 2x2 matrices. The symmetric matri- 
ces make a three dimensional space that is conveniently 
parameterized by cylindrical coordinates: 

S(t,p,0) =t + (p cos 6)a x + (p sin0) a z p>0. (8) 

In this parametrization det S — t 2 — p 2 and positive ma- 
trices correspond to the cone t > p. Since ia y is the 
generator of rotation in two dimensions A = —ia y dip. 
Eq. ijHJ can be readily solved for the connection dip whose 
differential is the curvature F. One finds 



dip 



1 



i + C 2 



■dS, F 
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(i + C 2 ) 1 



■d(C 2 )Ad9, (9) 



where Q = t/ p. F and dip are invariant under scaling of 
S, as they must be. 

The total curvature associated with the area preserv- 
ing deformations is then J t2 _ 2=1 F = \ J d6 = n. This 
implies that in any single closed cycle (i.e. one with- 
out self intersections), the angle of rotation is at most it. 
When interpreted as magnetic flux, tt corresponds to half 
a unit of quantum flux. 

A geometric understanding of why the total curvature 
of orientation preserving deformations (t > p in Eq. (JSJ) 
is 7r comes from considering also the orientation reversing 
deformations (t < p in Eq. (JSJ). When both are consid- 
ered £ takes values in [0, oo) and then F coincides with 
the curvature of the canonical line bundle (Berry's spin 
half) on the (stereographically projected) two sphere. 
The Chern number of the bundle is 1 and the total curva- 
ture 2tt. The orientation preserving matrices correspond 
to the northern hemisphere, £ > 1, and orientation re- 
versing to the southern hemisphere, £ < 1, |{|. 

The cost function must include some measure of dissi- 
pation. For, without dissipation energy is a function of 
the controls and no change in energy is associated with 
a closed loop. We consider two models of dissipation in 
isotropic media. Both lead to the same metric on the 
space of controls, namely 



(d£) 2 = Tr (S- 2 d(S 2 ) <g> S~ 2 d(S 2 )) 



(10) 



Consider a medium with viscosity tensor r\. The power 
due to dissipation is 



(11) 



where u is the strain tensor and u the strain rate. S is 
related to u by 2u = S 2 — 1. To see this recall that the 
strain is defined as the change in the distance between 
two neighboring points caused by a deformation. If the 
deformation is described by a metric g then g = 1 + 2u, 
where 1 is the metric associated with the undeformed ref- 
erence system . When considering linear deformation 
described by a symmetric matrix S, the resulting metric 
is g = S 2 (regarding the covariant components of g as 
the elements of a positive matrix). 

The space of (symmetric) 4-th rank isotropic tensors 
is two dimensional and spanned by the two viscosity co- 
efficients r\ and r)' 



V m = '///'.< 



v'g tJ g kl 



(12) 



We therefore find that in an isotropic medium 

P = 2 V Tr(g- 1 gg- 1 9 )+2r 1 'Tr (g^^Tr (g^g) (13) 

For volume preserving transformations the term multi- 
plying r/ vanishes and one is left with the first term alone. 
By choosing the unit of time appropriately one can take 
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i] = 1. This leads to the metric of Eq. I|10|) which is in- 
variant under congruence, g — > AgA*, for an arbitrary 
invertible matrix A. 

The dissipation in certain plastic materials can be rate 
independent. This is the continuum mechanics analog of 
the dissipation due to friction when one body slides on 
another |9(. In plastics, rate independence is a conse- 
quence of the Levy-Mise constitutive relation: s SX = Su, 
where s is the stress and 5X a scalar valued function ||. 
The constitutive relation is formally the same as for fluids 
[j| , and by isotropy, the dissipation must be a function of 
dt of Eq. i|10|) . If the material is memoryless, dissipation 
is additive with respect to concatenating paths and must 
be proportional to dl. 

Returning now to two dimensions, let us parameter- 
ize the area preserving transformation by: (t, p, 9) — 
(cosh(r/2),sinh(r/2),6'). The metric Eq. JTHJ gives 



(d£) 2 = {drf + sinh 2 r{ddf 



(14) 



This metric gives the hyperboloid the geometry of the 
pseudo-sphere i.e. it makes it into a surface of con- 
stant negative curvature —10. Geometrically, this cor- 



responds to embedding the hyperboloid t 
Minkowski space. 
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FIG. 1: Positive, area preserving linear transformation may 
be identified with the hyperboloid of revolution t 2 — p 2 = 1. 
Embedding the hyperboloid in Minkowski space gives it the 
structure of the Lobachevsky plane with constant negative 
curvature. 

The metric enables us to assign a scalar / to the curva- 
ture 2-form F of Eq. as the ratio between F and the 
area form of the pseudo-sphere, sinh r dr A dO. Similarly, 
a(r) is the scalar relating the length form sinh r d9 to the 
one-form dip. A calculation then gives 



/(r) 



1 



2 cosh r 



o(r) 



cosh r — 1 
sinh 2r 



(15) 



and is plotted in Fig. [21 The curvature is everywhere 
positive and it is concentrated near the origin, r = 0. It 
decays exponentially with t. This means that large de- 
formations are ineffective. We already know that small 
deformations are ineffective. This brings us to the opti- 
mization problem. 

The control problem is to find a closed path 7 in the 
space of deformations, starting at Sq = \/Q , (Qo is the 
initial quadrupole), which rotates the quadrupole by $ 
radians, with minimal dissipation. If the dissipation is 
rate dependent, one adds a constraint that the time of 
traversal is 1. 



FIG. 2: The scalar curvature, /(r), which can be interpreted 
as a magnetic field (left) and the 9 component of the vector 
potential a(r). / is exponentially localized near the origin of 
the control space (r, 9), while a has a maximum on the cir- 
cle coshr = (f>, (f> the golden ratio. a 2 /2 appears also as a 
potential in the effective 1-dimensional dynamics of the opti- 
mization problem. 



For viscous media 7 = (r(t), 9(t)) is then the solution 
of the variational problem 



1 /f 2 + 2 sinh 2 r 



Aa(r) sinh(r)6» ] dt = (16) 



where A is a Lagrange multiplier and 7(0) = 7(1) = 
So = S(to,6q). This is the (variation of the) action of 
a classical particle with charge A and unit mass moving 
on the hyperbolic plane (i.e. the pseudo-sphere) in the 
presence of a magnetic field /(r) given in Eq. (|15fl . 

Since motion in a magnetic field conserves kinetic en- 
ergy the particle moves at constant speed and the dissipa- 
tion, ^|7| 2 , depends only on the length of the path. The 
variational problem can therefore be cast in purely geo- 
metric terms: Find the shortest closed path starting at 
a given point which encloses a given amount of magnetic 
flux. The shortest path is evidently also the solution in 
the case that the dissipation is rate independent. 

Consider the family of isospectral deformations r = 
const which keep the the eigenvalues of S (or Q) constant 
while rotating its eigenvectors. We call these "stirring" , 
see Fig. |3J One can not stir an isotropic body, since its 
eigenvectors don't have well defined directions. 

Among the stirring cycles there is an optimal one 
which maximizes the rotation per unit length. From 
Eqs. (|9I14|) the flux to length ratio for stirring cycles 
is a(r) of Eq. 10. The function a takes its maximum 
at coshr s = <f>, t s w 1.061, see Fig. |3 Every cycle of 
the controls rotates by $ s — (2 — <f>)ir w .3827T radians, 
somewhat less than a quarter turn. 

To make full use of the optimal stirring cycle the ini- 
tial conditions must be right. This is the case for a 
quadrupole with Q = S(2t s ,6), arbitrary. With other 
initial conditions and for large angles of rotations, the 
optimal paths approach the optimal stirring cycle, linger 
near it, eventually returning to the initial point. This is 
shown in Fig. 0] for a half turn of an isotropic body. 

Since the magnetic field F(t) is rotationally invariant, 
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FIG. 3: A reference shape (left), assumed to have Q — 1, 
and four instances from the optimal stirring cycle (right). 
The configurations on the right are ordered clockwise with 
increasing 6, starting with the top left. Both the first and the 
last correspond to 9 = and are therefore related by a pure 
rotation. 
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FIG. 4: The optimal path in the space of shapes for a rotation 
by 7r of a body with initial quadrupole moment Q = 1. Polar 
coordinates (r, 6) are used to parametrize the plane. The 
path reaches exponentially close to the optimal stirring cycle 
at coshr = (f>, and winds around it twice before returning 
to the initial configuration. The optimal stirring cycle is not 
distinguishable from the envelope of the orbit in the scale of 
this figure. 



value of 9 when r = is meaningless. It then follows from 
Eqs. i|17[l and i|18fl that the optimal orbit, r(0), depends 
only on the ratio E/X 2 . Rescaling time properly, we can 
achieve A = 1 by relaxing the constraint that the time to 
complete a cycle should be 1. 

The key equation controlling the dynamics is Eq. I|18fl , 
which describes effective one-dimensional motion in the 
potential a 2 (r)/2. As shown above, a(r) has a maxi- 
mum at t s . Therefore, closed orbits which correspond 
to optimal paths have energy values E < a(r s ) 2 /2. The 
motion is quite simple: Trajectories leave the origin with 
a positive r, reach the turning point r t — a~ 1 (\/2E), and 
return symmetrically to the origin, thereby completing a 
cycle. The flux accumulated during a cycle is 

$ = [ a(r) sinh(r)M , (19) 
Jo 

where 6 is obtained from Eq. I|17l) . 

Increasingly longer orbits are obtained when E ap- 
proaches the separatrix energy a(r s ) 2 /2. The orbit cor- 
responding to the unstable equilibrium point at r = t s 
is the optimal stirring cycle, which is a J = orbit. The 
flux accumulated during a complete turn in parameter 
space, where 9 increases by 2ir, is hence bounded by <& s . 
Large values $ require many turns during which the or- 
bit approaches exponentially the optimal stirring cycle, 
see Fig. 0] A circular disc of playdough would therefore 
rotate by 7r, with no angular momentum, in about three 
cycles. 
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the angular momentum 

J = sinh r (o sinh r - Aa(r)) (17) 
is conserved. Conservation of energy gives 

0<f 2 =2E- (^r- +\cl(t)) , (18) 
\ sinh r / 

which reduces the problem to quadrature. The three con- 
ditions, t(1) = t(0), 0(1) = 0(0) mod 2vr, and the con- 
straint of enclosed flux determine the three parameters 
E, J and A. 

The initial condition r(0) = is special: By Eq. Ijl7|l 
the angular momentum J is forced to have the value 0. 
In turn, there is also one less condition to satisfy, as the 
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